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Abstract 

We study the t-J model in one dimension by numerically projecting the true 

ground state from a Luttinger liquid trial wave function. We find the model 

exhibits Luttinger liquid behavior for most of the phase diagram in which 

interaction strength and density are varied. However at small densities and 

high interaction strengths a new phase with a gap to spin excitations and 

enhanced superconducting correlations is found. We show this phase is a 

Luther-Emery liquid and study its correlation functions. 
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The t-J model was proposed to describe the dynamics of holes doped into a Mott insu- 
lating state |1]-§J. Even in one dimension, determining the complete phase diagram for this 
apparently simple model has proven to be quite formidable, and the ground state structure 
turns out to be far richer than initially suspected. In this paper we combine a variational 
approach with an exact ground state projection method to study the properties of this 
model. 

The Hamiltonian for the t-J model in one dimension can be written in the subspace of 
no doubly occupied sites as 

H = -tJ2i c la c i+i<7 + 4+i*Civ) + J H(Si ' S i+ i - jnin i+1 ). (1) 

ia i 

The model has been solved exactly only for J — > 0, where it is equivalent to the U — > oo 
Hubbard model, and J = 2t @||. In both cases the ground state at arbitrary density 
belongs to a broad class of interacting Fermi systems known as Luttinger liquids, which 
exhibit power law decay of correlation functions with exponents characterized by a single 
parameter Additionally, for very large J/t the attractive Heisenberg interaction term 

in (TJ) dominates the kinetic energy and the model phase separates. 

To obtain the rest of phase diagram of the t-J model, several numerical approaches have 
been used. For example, Ogata, et. al. |J have exactly diagonalized this Hamiltonian on 
a 16 site ring and find the model behaves as a Luttinger liquid for all values of J/t below 
a critical value at which phase separation occurs. They hypothesized that a third phase of 
bound singlet pairs may separate the other phases at very low density but were unable to 
resolve this phase with such small system sizes. 

In this paper we employ a Luttinger liquid variational wave function to approximate the 
ground state of the one dimensional t-J model [|I0|-|13|| and then use a numerical projection 



technique to extract the true ground state from this trial state. With these methods, we 
can accurately investigate much larger systems than attainable by previous techniques. We 
confirm that the t-J model has a Luttinger liquid ground state for most of its uniform density 
phase diagram, and in this region the ground state is well described by the trial state. At 



small densities, however, we find a third phase separating the Luttinger liquid and phase 
separated states. This phase behaves as a Luther-Emery liquid, exhibiting a gap to spin 
excitations and enhanced superconducting correlations. 

In previous work [0 we studied the ground state with a Luttinger liquid trial state 
written in the subspace of no doubly occupied sites as a Jastrow Slater wave function 

= 111^1" S(rj)S(rl) (2) 

i<j 

where S(ri) = Det[e tk ^ rt } is a Slater determinant of single particle plane wave states and 
dij = sin(7r(rj — rj)/L) for a system of size L. The Jastrow factor Yl^j \dij\ u in (||) modulates 
the wave function by the distance between all pairs of particles raised to the power u, taken 
as a single variational parameter. Positive values of v induce a smooth correlation hole 
between all particles, while negative values provide an attractive correlation competing with 
the Pauli repulsion. For v < — 1/2 this attraction overcomes the statistical repulsion, and 
phase separation occurs. The long range nature of this Jastrow factor generates the Luttinger 
liquid behavior of the wave function [0,0 . This wave function has been considered in two 



dimensions where it also exhibits an algebraic singularity at the Fermi surface [15 JIB | . 



Applying (^) to the t- J model one finds the optimum value of the variational parameter 
v varies continuously with interaction strength and density over most of the phase diagram 
prior to the critical J/t for phase separation. However, at small densities we found a third 
region separating the Luttinger liquid and phase separated states where the optimized pa- 
rameter is pinned at the critical state v = — 1/2. At this point the many body system in 
the trial subspace has infinite compressibility, which physically cannot extend for a range of 
interaction strengths. One concludes that the true ground state here likely lies far from our 
variational subspace. We would like a systematic way of both checking the accuracy of the 
trial state where we think it is doing well and determining the exact ground state in this 
third region. 

In this work, we start with the optimized trial state (Q) and project it onto the exact 
ground state numerically [p7H20|| . A series of increasingly accurate approximants to the 



ground state is generated by \p) — (H — W) p \^ v ) were H is the Hamiltonian and W is 
a numerical constant. These states approach the true ground state for large p provided 
\Eq — W\ > \Ei — W\ for all excited states E^. For the t- J model with J > we may choose 
W — 0. In principle this method can be used to project any trial state not orthogonal to the 
ground state, but good initial states give faster convergence and smaller statistical errors. 
To evaluate ground state expectation values of an arbitrary operator, we calculate 



{ ° p) (p\p) (¥„|J**|*/ (3) 

and take the large p limit. For sufficiently large powers, the scaling of fl3|) is dominated by 
the contribution from the first excited state overlapping the trial state. Thus we can write 

(E p ) = E + 5exp(-2Ap) + --- (4) 

with exp(— A) = \E\ — W\/\Eq — W\. An operator not commuting with the Hamiltonian 
has an additional cross term: 

(O p ) = O + 5 1 exp(-Ap) + 5 2 exp(-2Ap) + ■ • • (5) 

We use the convergence of the energy (f|) to fix A, and then use @ to determine the ground 
state values of the rest of the observables. 

Traditionally (§) has been calculated using a hybrid of two numerical techniques. First 
the trial wave function \l/ is sampled with Variational Monte Carlo to give an ensemble of 



configurations \a) Then for each |a) the product H p is sampled stochastically 



using a method similar to the Neumann Ulam matrix method [p3|| . The products are sand- 
wiched to evaluate (H p OH p ) or the normalization (H 2p ). This approach throws away much 
information, specifically the details of the intermediate states in the evaluation of each H p . 

We developed a much more efficient algorithm for evaluating (^) by combining the two 
operations. In usual Variational Monte Carlo a new configuration \j3) is chosen from a 
previous configuration |a) with probability 



P™ c = min 1, 



^ 2N 



(6) 



After many transitions, this leads to a distribution of configurations proportional to |^/«| 2 . 
If new configurations are instead chosen with the probability 

Pa-^p = —-tr H Pc* (7) 



z v& 



with 



= E W- H t»« ( 8 ) 
the distribution for a configuration |a) approaches l^/z^. This method of generating 
new configurations is the same used to evaluate the products H p , so the operations can be 
combined. 

The algorithm improves on the traditional approach in two ways. When evaluating a 
diagonal expectation value, such as (n(r)n(O)) or (S z (r)S z (0)) , our method evaluates a new 
(ty p\H p OH p \ i $> a) at every step of the random walk, so calculations of different powers p 
require the same amount of time. Additionally, for any expectation value, an arbitrary 
number of different values of p may be calculated in parallel. The only disadvantage of our 
approach is that ergodicity is violated as J — > 0, and the old method must be used in this 
limit. Both methods offer an improvement to Greens Function Monte Carlo in that exact 
correlation functions can be calculated pi ]. Since statistical errors grow with increasing p, 



we generally chose the maximum power to be 10 times the system size. 

The phase diagram of the t-J model determined by our projection technique is shown 
in Fig. [1]. We see that three distinct phase occur. For small J/t, the ground state is a 
Luttinger liquid with spin correlations dominating the long range behavior. Increasing J 
suppresses these correlations, and the ground state passes through the Fermi liquid point 
of the Luttinger liquid spectrum at the dashed line. Above this line the Luttinger liquid 
has dominant singlet pairing correlations, and for very large J/t the ground state is phase 
separated P, p!2|j25| , [26| . As will be reported in detail elsewhere, in the Luttinger liquid regime 
the trial state (0) approximates the exact ground state very well. 

In this work we see clear evidence of a new Luther-Emery liquid phase (labeled "Spin 
Gap") separating the Luttinger liquid and phase separated states at small densities 



Unlike all Luttinger states, this new phase exhibits short range spin correlations, and thus a 
gap to spin excitations, while both charge and singlet pair correlations decay algebraically. 
Physically one can view the Luther-Emery liquid as a translationally invariant coherent 
quantum fluid of bound singlet pairs. The pairs are correlated and can be treated at a 
simple level as an interacting fluid of hard core bosons. 

Luther-Emery states have been observed in diluted spin models that exhibit gaps in 
the saturated state, such as the t-J model with Ising anisotropy or the next nearest 
neighbor t-J model p9fl . Additionally this phase is present in the t-J- V model at quarter 



filling pOfl . This work provides the clearest evidence to date of the spontaneous formation 
of a Luther-Emery state by doping a gapless parent state. 

A sample spin correlation function in the Luther- Emery phase is plotted in Fig. (Q) with 
the correlation function obtained from the unprojected trial state shown for comparison. 
The variational function exhibits the linear behavior at small wave vectors characteristic 
of Luttinger liquids, while the exact function is quadratic at small k and analytic at all 
wave vectors, consistent with exponentially decaying spatial correlations. We calculate the 
boundary between the Luttinger and Luther-Emery states by the crossover from linear to 
quadratic behavior at small wave vectors. 

More definitive evidence of Luther-Emery behavior can be seen in the superconducting 
correlation function plotted in Fig. |3|. The exponents of the correlations functions in both 
Luttinger and Luther-Emery liquids that decay with power laws can be characterized by a 
single parameter K p > pl| , |32| . The non-oscillatory part of singlet pair correlation function 



decays as 

(b\r)b{0)) ocr- x (9) 

where b{r) = ^(c r |C r+1 j — c r jC r+ if). For Luttinger liquids = l+K' 1 while Luther-Emery 
liquids have Xl-e = 

In Fig. |3], b(k) diverges logarithmically with system size as k — > in our trial wave 
function, which represents the strongest divergence possible in a Luttinger liquid state. 



However, the true ground state in this region apparently exhibits a much stronger cusp |T9 . 

Using a finite size scaling analysis of the divergence in b(k — > 0) in the projected state, 
we can determine the value of this exponent A ||25|| . A plot of A showing the transition from 
Luttinger to Luther-Emery liquid behavior at density 5 = 1/6 is shown in Fig. [|. In the 
Luttinger regime, A is bounded from below by 1, but this bound is clearly violated as the 
Luther Emery state is entered. A continuous variation of A with J as found in this data 
would imply a discontinuous jump in K p . 

It is interesting to note that non-interacting hard core bosons have A = 1/2, so our singlet 
pairs have residual repulsive interactions for J < 2.65 in the Luther-Emery state, while at 
higher J the hard core nature of the pairs competes with an effective attractive interaction 
| 33fl . The attraction from the Heisenberg term in ([]]) in this regime is strong enough to bind 



singlet pairs but still insufficient to cause macroscopic phase separation. 

Chen and Lee proposed a variational state for this region by Gutzwiller projecting a sea 
of non-interacting bound singlet pairs . Their wave function corresponds to a K p = oo 
Luther-Emery state, the critical point of the verge of phase separation which exhibits a 
macroscopic superfluid density. Their calculations of the boundaries of the spin gapped 
regime agree remarkably well with ours except at the boundary to phase separation, which 
they find occurs at higher J/t. One may speculate that a potentially more accurate trial 
state could be generated by correlating the pairs with a Jastrow factor similar to (0). This 
state would exhibit generalized Luther-Emery behavior with arbitrary K p . 

In summary, we have investigated the ground state properties of the t-J model in one 
dimension using a numerical technique to project the exact ground state from a variational 
Luttinger Liquid trial state. We find the model has a surprisingly rich phase diagram. At 
lower interaction strengths the variational wave function accurately describes the Luttinger 
liquid phase, and at very large J/t the model phase separates. However, one finds these 
phases are separated at low density by a Luther-Emery quantum dimer liquid phase with 
short range spin correlations and enhanced superconducting correlations. 
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FIGURES 

FIG. 1. The phase diagram of the t- J model as determined in this paper. The Luther-Emery 
state is the region labeled "Spin Gap". The dashed line indicates K p = 1, the Fermi liquid phase. 
Below this line the Luttinger liquid has dominant antiferromagnetic correlations, while above this 
line singlet pair correlations decay with the smallest exponent. The phase separation boundary is 
determined by the divergence of n(k — > 0), the lower Luther-Emery boundary by the behavior of 
S(k — > 0), and the Fermi liquid line by S(k — > 2kp). All systems contained at least 100 sites and 
10 electrons and holes, so phase boundaries cannot extent to the extreme densities. The dotted 
lines are extrapolations. 

FIG. 2. The spin-spin correlation function for J/t = 2.8 and density n = |. The optimized 
variational wave function has linear behavior at small wave vectors while the exact spin correlation 
turns on quadratically in k. The system contains 20 electrons on a 120 site lattice. 

FIG. 3. The singlet pair correlation function at J/t = 2.8 and density n = |. The k = cusp 
is greatly enhanced in the exact ground state. The system contains 10 electrons on 60 sites. 

FIG. 4. The scaling of the exponent of the k = superconducting cusp with interaction 
strength at density n = 1. The transition from Luttinger to Luther-Emery liquid states occurs 
at J/t « 2.3 and the system phase separates at J/t « 2.9. Luttinger liquids require A > 1, and 
non- interacting hard core bosons have A = ^. 
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